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Abstract. If cold dark matter elenientary particles form a Bose-Einstein 
condensate, their superfluidity may distinguish them from other forms of cold 
dark matter, including creation of quantum vortices. We demonstrate here that 
such vortices are favoured in strongly-coupled condensates, while this is not the 
case for axions, which are generally presumed to form a Bose-Einstein condensate 
but are effectively non-interacting. 



Introduction 

Suggestions have appeared in the literature that cold dark matter (CDM) may 
be in the form of a Bose-Einstein condensate (BEC), including axionic and other 
forms of CDM. This has important implications for the physics of structure for- 
mation, notably at small scales where one expects significant deviations from 
the more standard CDM, due to the superfluidity exhibited by BECs. In fact, a 
prime motivation for considering self-interacting BECs for CD M was their ability 
to p roduce galactic halos with constant density cores, see e.g. iGoodmanI (|200m 



and lPeeblesI (I2OOOI '). The corresponding profiles may then bet ter agree with ob- 



served rotation curves of dwarf an d LSB galaxies, see e.g. iBohmer k, Harkj 



(l2007l ). iBohmer. Martins fc Saluccil ^M)- The inner structure of halos may 



also depend on the angular momentum distribution of the CDM particles. Lab- 
oratory BECs are known to develop v ortices when rotated with a sufficient an- 
gular velocity. iDuffv &: Sikivid ( 20081 ) report that certain fine-structure in the 



observed inner mass distribution of the Milky Way can be explained only if 
the infalling dark matter particles from which such galactic halos formed had 
a net overall rotation, causing a 'tricusp' caustic ring of dark matter in that 
case. For standard, no n- interacting CDM mo dels, however, one expects infall 
to be irrotational, while Sikivie Sz Yang ( 20091 ) argue that axionic dark matter. 



as a BEC, may form vortices leading to net overall rotation. As such, they 
suggest. Milky- Way observations may already have detected the signature of 
axionic CDM. Vortices have also b een postulated for s t rongl y-coupled BECs 
involving ultralight scalar particles (j Silverman &: Mallet ( 20021 )). The question 



of whether an angular velocity sufficient to create vortices occurs in BEC dark 
matter cosmologies has not yet been answered, however. We address this point 
here by calculating the critical angular velocity for vortex creation in a simple 
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model of BEC/CDM galactic halos and comparing the result with the angular 
velocity expected from cosmological N-body simulations of CDM. 



The Model 

We shall use an energy argument to derive the critical angular velocity for vortex 
creation in a rotating, self-gravitating BEC halo by finding the angular velocity 
above which the energy is lowered by the presence of a vortex. For the un- 
perturbed equilibrium state, we model the BEC halo as an oblate Maclaurin 
spheroid - a homogeneous ellipsoid of mass density p = mn, with semi-axes 
(a, 6, c) along {x,y,z) such that a = 6 > c, of total volume V = Aira'^c/S, 
uniformly rotating with angular velocity Q = (72, with gravitational potential 
<I>(r, z) = TrGp (j4i(e)r^ + A3(e)z^) in cylindrical coordinates {r,z). The func- 
tions Ai{e), A^i^e) depend on the excentricity e, defined via e = y^l — {c/a)"^, 
and there is a family of solutions par a.meterized by 0/0(7 , or, e quivalently, e, 
related to each other according to (see lBinnev &: Tremaind ( 19871 )) 



= 2[^i(e)-(l-eV3(e)], (1) 



where Oc = yJnGp is a characteristic gravitational angular frequency. For 
O = 0, e = 0, while e must not exceed 0.9529 (ibid.). 

We describe these self-gravitating BEC halos of ellipsoidal shape with vary- 
ing degrees of rotational support by self-consistently coupling the Gross-Pitaevskii 
(CP) equation of motion for the complex scalar wavefunction ^{r, t) to the Pois- 
son equation, where |'(/'P(r,t) = n, the number density of particles of mass m: 

dih 

= AV' + lm^ + c/l^/;!^ -/x)V', A$ = 47rGm|^/;|^ (2) 

We assume that halos are comprised of N particles in the condensed state de- 
scribed by ^, and so f^^ = N, which determines the CP chemical potential 
//. BEC/CDM, like standard CDM, is assumed to interact so weakly with other 
matter and radiation, once its abundance is fixed in the early universe, that 
we can neglect all other, non-gravitational couplings. However, BEC dark mat- 
ter is self-interacting, as described by an effective interaction potential g|V'|^/2 
with coupling constant g = Airh'^as/m, where is the 2-body scattering length. 
[Since BECs with negative scattering length are not stable in the context de- 
scribed here, we consider only g > 0.] 

In a frame rotating with velocity Ct, our system is stationary, so the deriva- 
tive with respect to time in equ.dJ]) vanishes. The equation of motion in this 
frame is then given by equ.([2]) with an additional operator Vrot = ■ L on the 
right-hand-side, where L = —ihr x V, and it is also understood that the respec- 
tive variables and quantities are in the new frame. Incompressible Maclaurin 
spheroids are approximate solutions of this system of equations Q Such station- 



^ Uniform density is a more realistic approximation for BEC halos than for the cuspy halos 
of standard CDM; equ. © - with or without rotation - favours a flat core over the cuspy 
r~^-profile found by N-body simulations of standard CDM. 
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ary systems can then be studied via the corresponding GP energy functional, 
given by 



|!.|v^|2 + |ci,|v,|2 + + mrn ■ (r X w) 



(3) 



We shall use this equation to determine at which angular velocities the 
presence of a vortex is energetically favoured. To this aim, we decompose the 
wave function ^ into a vortex-free part and a part which carries the vorticity, 
thus splitting the energy functional in equ.([3]) so as to compare the respective 
energy contributions, with or without vortices, more easily. To derive an analyt- 
ical result with as much generality as possible, we consider the following ansatz. 
The vorticity-part is a d-quantized straight vortex, modeled as a funnel-shaped 
tube along the rotation- axis, with core radius s, where the density outside of 
the core is given by the unperturbed profile of the Maclaurin spheroid, whereas 
the density in the vortex core region drops to zero at the center. 



Results 

When our wave function ansatz is inserted in equ.dS]), there is a critical angular 
velocity Jlc above which the energy is lowered by the presence of a vortex. Since 
our ansatz leads to an energy greater than or equal to that which would result 
if the 'real' wave function were used in place of the ansatz, this means vortex 
creation is energetically favoured, in general, if > Qr. The general expression 
for Op (for brevity not shown here; for derivation see iRindler-Daller &: Shapirol 



([20101)) depends on the core radius s, which is of the order of the healing length 
defined by equating the energy contributions from the kinetic term and self- 
interaction, so = /{2pg). In the strongly-coupled regime, as g — )• cxd, ^ — )• 
and s — 7- 0, respectively, and the critical angular velocity has a leading-order 
term which diverges logarithmically, Oc — ^QAidln.^, just as for laboratory 



condensates (see e.g. iLundh et al.1 (|l997l ^). Here we define ^qm = h/{ma'^), the 



characteristic angular frequency such that every particle contributes an amount 
h to the total angular momentum of the uniformly rotating Maclaurin spheroid 
(|L| = Nh). On the other hand, if s or ^ are comparable to the system size 
a, this logarithm is subleading. We note that Qc is a monotonically increasing 
function of the winding number d, so we restrict our consideration to the lowest 
Qc, for which d = 1. Then, the co re radius can be replaced by the healing length 
(see iPitaevskii &: Stringarl ( 20031 )). and the critical angular velocity (in units of 
Oc) becomes 



no nc ' 12 ' 2 yriQMj \ Q \U ' 9 

(4) 

We determine whether a given set of BEG parameters (m, g) makes the 
spheroid rotation velocity equal the critical value above which a vortex forms, 
O = Qc, by combining equs.([l]) and (jH, for a given e. We observe that Qqm/^g = 
niH/m anda/^ = {g/gminY^'^, where m// = h/[a?{'KGpY/'^] and^mm = f? I{2pc? 
the latter of which depend only on halo parameters. The solution for a given e 
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is a curve in the {m/niH , g/ gmin)-P^Sine, for which there is a minimum allowed 
value of g/gmin ^ 1- For each curve, no vortex is allowed for parameters in the 
space above the curve. According to equ.dH), goes to infinity as g' — >• 0. This 
is the case for axion dark matter, for which the coupling is so weak as to be 
effectively zero. 

We determine the e- values of interest for BEC/CDM halos as follows. The 
superfluidity effects of BEC dark matter which distinguish it dynamically from 
standard CDM are mostly limited to the internal structure of our halos, while 
larger-scale structure formation is otherwise the same. The latter is responsi- 
ble for the tidal torques that give a halo its angular momentum. Cosmological 
N-body simulations of the CDM universe show that halos form with a net an- 
gular momentum such that the dimensionless ratio A = L\E\^/'^ /GM^/"^ ^ which 
expresses their degree of rotation al support, has value s in th e range [0.01,0.1] 
with median value 0.05 (see e.g. iBarnes fc Efstathioul ( 19871 )). where L is the 



angular momentum, E is the halo binding energy, and M is the total mass. For 
our spheroids, corresponds roughly to the ratio of rotational kinetic energy 
to gravitational potential energy, and can be expressed in terms of e only, 

1/2 



6 arcsin(e) / 3 , 3\/l — 



5y^ e \ 2e2 2earcsin(e) 



In what follows, we take three representative values for A, (0.01, 0.05, 0.1), which 
correspond to e = (0.051, 0.249, 0.464), respectively. For given values of the halo 
mass density and eccentricity, the angular velocity is fixed. For example, a 
Milky- Way-sized halo with M = 10^^ Mq and a radius of -R = 100 kpc, where 
the corresponding Maclaurin spheroid with the same mean mass density has a 
semi-axis given hy R = a(l — e^)^/^, has a density of about 10~^^ g/cw? and 
n ~ 10"^^ rad/s for A = 0.05, so ~ 0.18. For such a halo tuh ~ 10"^^ g 

and gmin ~ 10"^® erg cm^. In Fig.l, we plot m/niH versus g/gmin for the above 
A-values, independent of halo size. 

For BEC halos, the condition of gravitational equilibrium restricts the val- 
ues of m and g to another curve in the {m/mH,g/gmin)-'9^8ja.e. In the strongly- 
coupled regime, non-rotating BEC halos are just (n = l)-polytropes, for which 
the size R is related to the B EC parameters according to R = T:[f\^as/ {Gw?)Y/'^ 
(see Bohmer &: Harkol ( 20071 )). which we translate into our language as 



J 8(l-e2)i/3 



(5) 



The resulting relationship between m/niH and g/gmin is plotted in Fig.l, as 
well. For CDM halo A-values, the degree of rotational support is small enough 
that equ.dS]) should still be a good approximation. There is almost no sensitivity 
to e, so the respective curves for different e lie on top of each other. Since the 
BEC parameters which satisfy ([5]) are all below our critical curves for which 
= r^c, BEC/CDM halos, in general, will typically form vortices. 
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Figure 1. Dimensionless BEC particle mass m/niH vs. coupling strength 

g/ grain- Critical lines — ^Ic ioi \ — 0.01 (e — 0.051) (grey-dotted), 

A = 0.05 (e = 0.249) (black-solid), A = 0.1 (e = 0.464) (light grey-dotted); 

BEC halo (n = l)-poly tropes: lower- most curves (grey-solid) for the same 

e-values 



Conclusions 

While axions, with an effectively zero coupling, apparently do not form vortices, 
vortices will be created, in general, for BEC/CDM halos in the strongly-coupled 
regime. As such, previous BEC models of halo mass profiles, which do not ac- 
count for their presence, should be revised, especially at the centers where vortex 
creation causes the density to drop. 
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